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Two-loop QED corrections to the decay rate of para- 
positronium (p-Ps) into two photons are presented. We 
find r(p-Ps 77) = 7989.64(2) /is"\ The non-logarithmic 
O(a^) corrections turn out to be small, contrary to some ear- 
lier estimates. We speculate that the so-called "orthopositro- 
nium lifetime puzzle" will not likely be solved by large QED 
corrections. 
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Positronium (Ps) , the simplest known atom, is an ideal 
system to test quantum electrodynamics (QED) of bound 
states. The spectrum and lifetimes of Ps states are, at 
least in principle, calculable within QED with very high 
accuracy. Hadronic effects, which in other atoms limit 
the attainable theoretical precision, are suppressed by 
the small ratio of electron and hadron masses. 

The lifetimes of the triplet and singlet ground states 
(respectively orthopositronium and parapositronium) 
have been subjected to very precise theoretical and exper- 
imental studies. Theoretical predictions for parapositro- 
nium (p-Ps) and orthopositronium (o-Ps) decay rates 
into 2 and 3 photons respectively, can be expressed as 
expansions in the fine structure constant a: 
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are the lowest order decay widths of the p-Ps and o-Ps, 
respectively, and the ellipses in eqs. (|l],^ denote unknown 
higher order terms which we will neglect in our analysis. 
Corrections of 0{a) were calculated in ^ for p-Ps. For 
o-Ps the most accurate result was obtained in where 
references to earlier works can be found. The logarithmic 
two-loop correction was found in Q for o-Ps and in 0| 
for p-Ps. The leading logarithmic correction at three 
loops was computed in Some partial results on the 
0(0!^) corrections for both p-Ps and o-Ps can be found 
in but complete values of Bp o have not been 

obtained so far. 

Using eqs. (|l|,H,H), one obtains the following theoretical 
predictions for the lifetimes 



r*'^^/ = 7989.42 Ais-i + r(,"'Sp (^) ' 

^.theory ^ 7,0382 ^iS"! + T^B,, (^)' 
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How do these predictions compare with experiments? 
For p-Ps the most recent result |10| , 



rj^'lPp^ = 7990.9(1.7) ^JLS-\ 



(6) 



is in good agreement with (|j) if Bp is not too large. 

For orthopositronium the situation is not clear. Precise 
experiments of the Ann Arbor group |TT|Jl^ ] found 



r^1Ppg(gas measurement) — 7.0514(14) /is ^, 
r^^'^Pg (vacuum measurement) = 7.0482(16) /is~^, (7) 

which, for Bo = 0, differ from (H) by 9.4(t and 6.3cr re- 
spectively. This apparent disagreement of experiment 
with theory has been known as the "orthopositronium 
lifetime puzzle." It should, however, be noted, that a 
more recent Tokyo result |13| , 

r°'lPp^(Si02 measurement) = 7.0398(29) /iS-\ (8) 

agrees with the theory prediction if Bq is not too large. 
Both Tokyo and Ann Arbor groups are working to im- 
prove their results. 
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Should future experimental efforts confirm the Ann Ar- 
bor results (0), in disagreement with the QED prediction 

, the orthopositronium lifetime puzzle could be solved 
if Bo turns out to be unusually large, e.g. ^ 250 for the 
vacuum measurement. Alternatively, one might specu- 
late that some "New Physics" effects such as o-Ps decays 
involving axions, millicharged particles, etc., cause the 
excess of the measured decay rate over the QED pre- 
dictions. Some of those exotic scenarios seem to have 
already been excluded by dedicated experimental stud- 
ies. (For a review and references to original papers see 
e.g. 0.) 

It was anticipated that a full two-loop calculation 
might first be done for p-Ps. In fact the relative theoret- 
ical simplicity of p-Ps motivated the most recent lifetime 
measurement pO| . In this paper we present a complete 
calculation of Bp. Our result permits the rigorous test 
of bound state QED envisioned in . We find that the 
coefficient of the two-loop non-logarithmic term in (Q) is 



Bp = 5.1(3), 



(9) 



and the theoretical prediction for the p-Ps lifetime be- 
comes 



ptheory ^ 7939 64(2) ^s"!. 
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Below we briefly discuss some details of our calculation. 
The decay width of p — Ps ^ 27 can be written as 
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where P is the four-momentum of the p-Ps, 4>{p) is its 
wave function, and A{X, p) is the annihilation amplitude 
of an e'^e^ pair into a pair of photons with polarization 
{A}. 

In the non-covariant perturbation theory the on-shell 
amplitude of the process e"^e^ — s- 77 reads 
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Here Ep = ^ m? + p^ ; p and fc are electron and photon 
three- momenta in the p-Ps rest frame; and A±(p) are the 
projectors on the positive and negative energy states. 
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To leading order one can neglect the small momenta p ^ 
ma compared to m and |fc| ~ m. One finds the following 
leading order amplitude 



Alo = ^(de2)(afc)(aei). 
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and the lowest order decay width 
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Higher order corrections to Flo will be calculated us- 
ing non-relativistic QED (NRQED) ||l5| with dimensional 
regularization |16|. We divide up the corrections into 
three parts: 
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where Bp^^ 



^'"'^ is the contribution of the one-loop ampli- 
tude squared and B^^''^- are the hard and soft con- 
tributions. The hard corrections arise as contributions 
of virtual photon momenta fc ^ to. Their effects are de- 
scribed by adding operators containing 6(r) to the non- 
relativistic Hamiltonian. The technical challenge is to 
compute the Wilson coefficients of those operators. For- 
tunately those coefficients can be obtained using any con- 
venient external states. In particular one can compute 
them for the electron and positron at rest. It is impor- 
tant to employ dimensional regularization, so that one 
avoids the necessity of subtracting the Coulomb singular- 
ities from box graphs. On the other hand, the soft contri- 
butions come from the region of virtual photon momenta 
of the order of A; ~ ma and are sensitive to bound state 
dynamics. The actual calculation is briefly described be- 
low. 

The square of the one-loop amplitude is easily obtained 
from the one- loop result: 
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Bp^''^ consists of three types of contributions: vacuum 
polarization insertions in the photon propagators, light- 
by- light scattering diagrams, and two-photon corrections 
to the annihilation amplitude. 



B 
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^hard(yp) ^ 5hard(LL) + ^"^^(77). (18) 



Vacuum polarization insertions into the one-loop graphs 
(an example is shown in Fig. |l|(a)) were computed in 

r7i,ra. 



^hard(yp) = 0.4468(3). 



(19) 



Light-by-light scattering contributions (for examples see 
Fig. ^b,c)) are more difficult to compute because of their 
imaginary parts which make numerical integration un- 
stable. We computed them analytically, by formally as- 
signing a large mass M to the internal fermions and ex- 
panding in a; = m/M. The resulting series converge so 
well that several terms are sufficient to find the result at 
x=l: 



B^'"'^{IA.) = 1.28(13). 



(20) 
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The most difficult class of effects comes from the two- 
photon corrections, examples of which are shown in 
Fig. I^(d-f). We proceed in the following way: com- 
bine propagators using Feynman parameters; perform 
momentum integrations analytically; extract ultraviolet 
(UV) and infrared (IR) divergent pieces; integrate nu- 
merically over typically 5 (in some cases 6) Feynman pa- 
rameters in the finite expressions. 

Extraction of UV divergences is relatively simple - they 
show up, roughly speaking, as singularities in the overall 
factors rather than as divergent integrals over Feynman 
parameters. 

Calculation of IR divergent diagrams is more demand- 
ing. For each such diagram we add and subtract a simpler 
diagram so that the sum is IR finite, and the subtraction 
can be calculated analytically. As an example let us con- 
sider the diagram shown in Fig. |l|(d) . The IR singularity 
in this diagram appears when momenta of the virtual 
photons are small. To suppress the contribution of this 
region we rewrite the propagator of the t-channel electron 
by subtracting its value and a suitable number of deriva- 
tives when both virtual photon momenta are zero. This 
difference leads to an infrared finite expression which can 
be computed numerically. To compensate this subtrac- 
tion we must add the same diagram with the t-channel 
propagator replaced by a constant, l/(2m^). Such dia- 
grams were studied in [ p^ and can be computed analyt- 
ically. 

Similar tricks are used to compute all other IR di- 
vergent diagrams, although the subtraction procedure is 
more tedious in case of diagrams with stronger singular- 
ities, like the planar box in Fig. ||(e). 

Adding all two-photon diagrams and summing up nu- 
merical errors of individual diagrams in quadrature we 
find 



Bf{AA) = 
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The logarithm of the dimensionfuU parameter m arises 
from the expansion of the overall factor m~^' and van- 
ishes in all physically meaningfuU ex pre ssions. 

The sum of eqs. (|l|), and (|2l|) gives the total 

hard correction 



^hard 
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27r2 In m- 40.46(30). 
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To calculate the soft scale contributions, one should 
account for the relativistic corrections to the annihilation 
amplitude (A A) e+e^ — s- 77 and to the positronium wave 
function (WF): 



= Bf\AA) + Bf\WF). 



(23) 



For the annihilation amplitude correction, one expands 
the on-shell amplitude ( p^ ) to relative order Oip^ /m?). 
Although the resulting integral is linearly divergent, us- 
ing dimensional regularization one finds a finite result 
(see [p3 for a discussion of this effect): 



(24) 



Relativistic corrections to the positronium wave func- 
tion can be computed using the Breit Hamiltonian. Since 
we regularize all divergences dimensionally, we need the 
Breit Hamiltonian in d-dimensions derived in ||20| . Its 
projection on the S'-states can be found in eq. (39) of 
that paper. Performing calculations similar to those de- 
scribed after that equation, we find the wave function 
correction to the decay rate: 



B^°^\WF) H- 27r2 In - = — -h 2tt'^ In — 
^ a 2e 
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where on the LHS we have separated the logarithm, to 
be consistent with the division of corrections introduced 
in eq. §). 

The sum of the corrections to the annihilation ampli- 
tude (|2j) and to the wave function (|2^) gives the final 
result for the soft contributions, 
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We note that this partial result cannot be directly com- 
ared to the soft corrections found in a previous study 
since a different regularization scheme was employed 
there. 

The final result in eq. (^), Bp — 5.1(3), is obtained as 
a sum of the square of the one- loop corrections (|l^) , and 
the genuine two-loop hard ( p^ ) and soft (|2^) contribu- 
tions. We can now present a theoretical prediction for 
the 2-photon width of parapositronium with the 2-loop 
accuracy: 
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Our final result ( p7| ) agrees well within 1 a with the most 
recent experimental result, eq. (||). 

p-Ps can also decay into 4 or more photons. Those 
effects increase the p-Ps width by about 0.01 fis^^ (see 
pl| and references therein). 

The coefficient of the non-logarithmic (a/ir)'^ term in 
eq. ( PtI) is rather small, due to an almost complete can- 
celation between the soft and hard corrections. As we 
have already mentioned, only the sum of the two is regu- 
larization scheme independent and hence unambiguous. 
For this reason, it is likely that the cancellation between 
soft and hard pieces is not accidental. Scheme and gauge 
independent corrections, for example eqs. 




seem to indicate that "natural scale" of the 0{a 
rections is [several units] x (q/tt)^. 

The result of our calculation. Bp = 5.1(3), is much 
smaller than the estimate 40 ± 20 given in Q]. This 
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discrepancy can be traced back to the discussion after 
eq. (26) in the first paper of [Q. It seems that the impact 
of short-distance (hard) corrections was underestimated 
there. As the division of contributions into soft and hard 
pieces is regularization scheme dependent, it is clearly 
dangerous to draw conclusions about the complete result 
on the basis of only one of those parts. 

Having for the first time a complete two-loop correc- 
tion to a QED bound state decay, one is tempted to 
speculate about the size of such corrections to the or- 
thopositronium lifetime. Although nothing can be said 
rigorously, one could argue that most known second order 
effects have similar order of magnitude for o-Ps and p- 
Ps. A possible enhancement may be due to the larger (by 
about a factor of 3) number of Feynman diagrams con- 
tributing to o-Ps decay compared to p-Ps decay (this is 
already seen in the magnitudes of the 0{a) corrections). 
Unless this enhancement is dramatic for the complete 
0{a^) corrections to the decay rates, the theoretical pre- 
diction for the o-Ps lifetime will remain distant from the 
experimental results in eq. (0). It is therefore extremely 
important that the 3-photon decay of o-Ps be further 
studied experimentally. 
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FIG. 1. Examples of two-loop hard corrections to p-Ps de- 
cay into two photons. 
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